Hysteretic creep of elastic manifolds 



Stefan Scheidl^ and Valerii Vmokur^ 
(a) Institut fur Theoretische Physik, Universitat zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 
(b) Material Science Division, Argonne National Laboratory, Argonne, IL 60439 

(2 August 1996) 

We study the dynamic response of driven systems in the presence of quenched disorder. A simple 
heuristic model for hysteretic creep of elastic manifolds is proposed and evaluated numerically. It pro- 
vides a qualitative explanation of the phenomenology observed in experiments on high-temperature 
superconductors. 
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An intense attention is attracted currently by the 
driven dynamics of vortices in high temperature super- 
conductors. This interest is motivated not only by the 
technological request, but also because vortex dynamics 
exhibits features which are generic for a diversity of phys- 
ical systems including charge density waves, Wigner crys- 
tals, dislocation motion, polymers, and interfaces driven 
in quenched disordered media. In addition vortices are 
a unique experimental object with tunable parameters. 
One of the fundamental features of driven dynamics in 
quenched disorder is the existence of the so-called de- 
pinning transition: at zero temperature the system can- 
not move below a certain critical force whereas it moves 
above. At finite temperatures the transition is smeared 
out but the depinning threshold still marks the crossover 
from strongly non-linear thermally activated dynamics 
dominated by pinning below the critical force to the vis- 
cous flow above M. 

Pinning of the vortex lattice is governed by an inter- 
play between thermodynamic and kinetic effects. One 
of the manifestations of such an interplay the peak ef- 
fect in the critical current, has been recently observed in 
YBCO crystals prior to vortex lattice melting Fur- 
ther investigations of the transport just below the peak 
effect temperature range revealed another striking fea- 
ture of the depinning transition in the vortex lattice: 
dynamic hysteresis in the I-V characteristic and aging 
effects ||. It was found that the disappearance of the 
vortex motion (pinning) upon decreasing the driving cur- 
rent and the depinning upon keeping the vortex system 
at rest and subsequent current ramping occurs at differ- 
ent threshold currents: depinning threshold appeared to 
exceed noticeably the pinning current of the descending 
I-V curve. Moreover, at an ascending depinning thresh- 
old jth the voltage jumps abruptly from zero to the finite 
value V(jth) of the descending branch of the I-V curve 
showing thus a strongly asymmetric hysteretic behavior. 
The value of the ascending threshold current was found 
to depend on the waiting time i w the system has spent 
at rest: jth grows with the t w increasing. 

The observed abrupt switching of the vortex motion re- 
sembles the switching hysteretic behavior in charge den- 



sity waves transport Q and can be attributed to the for- 
mation of the bi-stable dynamic state in the vortex mo- 
tion, i. e. the dynamic state where both weakly pinned 
and strongly pinned vortex motions can coexist. The 
formation of the bi-stable state is usually related cither 
to the possibility of plastic deformations (phase slips) of 
the inhomogeneously pinned soft vortex lattice and/or to 
the formation multivalued equilibrium state for pinned 
vortices. It has been recently demonstrated that plastic 
deformations can become a determining component of 
the depinning of the vortex lattice Q . However the ex- 
isting descriptions of plastic deformation are lacking in 
aging and memory effects which are most essential and 
characteristic for glassy dynamics. In this Letter we will 
concentrate on this glassy aspect of the depinning process 
and discuss plastic models elsewhere. 

We consider a purely elastic threshold depinning pro- 
cess modeling the vortex lattice as an elastic manifold 
and propose a phenomcnological equation describing a 
temporal evolution of the typical barriers controlling the 
creep motion. We demonstrate the existence of aging 
effects which give rise to the abrupt increase in the man- 
ifold velocity within the creep regime and to a dynamic 
hysteresis in the depinning process. We find that the 
threshold force grows with the increasing waiting time 
t w the system spends at zero force. Our results pro- 
vide a solid basis for a description of the experimentally 
observed time-dependent hysteresis of the vortex lattice 
threshold dynamic behavior || in terms of the aging ef- 
fect and predict the possibility of similar effects in all the 
diversity of dynamic systems that can be modeled as an 
elastic manifold in a random quenched environment. 

We consider an overdamped dynamics of the elastic 
manifold subject to the quenched disorder. The charac- 
ter of motion follows from the interplay between the (i) 
elastic forces that tend to keep the manifold flat, (ii) pin- 
ning forces arising from the disorder potential which pull 
the manifold into potential minima forcing thus rough- 
ening of the manifold, (iii) thermal forces of a heat bath 
at temperature T, and (iv) a driving force F which is 
homogeneous in space but which may vary in time. 

The structure of a manifold in equilibrium is deter- 
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mined by the balance of elastic and pinning energies and 
can be characterized by the difference correlation (rough- 
ness) 



C(L) = (lu(L)-u(0)}2)*e (L/L ) 2 t 



(1) 



of deformations u from a flat configuration. £ denotes 
the transverse correlation length of the disorder and Lq 
is the Larkin length. Disorder leads to a rough manifold 
with an exponent £ > below four dimensions Q. The 
manifold behaves as in a landscape of energy barriers of 
all sizes L with heights U (L) which scale like 



U{L) = U (L/L ) 



(2) 



The exponent 8 — 2(£ — £o) is related to the increase of 
the roughness exponent £ due to the presence of disorder 
compared to the thermal roughness exponent (o = (2 — 
d)/2. The height of the smallest existing barriers on the 
Larkin scale is denoted by Uq. 

The typical time to overcome a barrier of height U by 
thermal activation is 



T{U)=t eMU/T] 



(3) 



with a microscopic temperature dependent time scale t$. 
The presence of barriers with unbounded heights gives 
rise to a glassy behavior of the manifold with diverging 
relaxation times. The central point of our analysis is a 
self-consistent treatment of the the dynamics of the man- 
ifold: we question the rate of the change of the manifold 
position and configuration given by its structure C(L), 
taking into account that the barriers governing the dy- 
namics are determined themselves by the instant mani- 
fold configuration. 

Relaxation into equilibrium. We start with the relax- 
ation of the perfectly flat manifold, C(L) = 0, into the 
equilibrium in the absence of a driving force. After a 
time t the typical depth of the metastable state the man- 
ifold falls in (or, equivalently, the height of the barrier 
separating the typical metastable states) is 

U(t)=Tln(t/t ) (4) 

Therefore the manifold is equilibrated up to length scales 

LJJJ) = L (U/U ) 1/e (5) 

according to Eq. (ph. In other words the manifold is 
rough up to this length scale, but it is still smooth above: 



C(L) 



${L/L } 2 < L<L C 
£ 2 [L C /L ] 2C L>L C 



(6) 



Here the time dependence enters through L c = L c (U(t)) 
via Eqs. (@) and (g). In the course of time the manifold 
becomes rougher, the barriers increase and the dynamics 
becomes even slower. 



Stationary driven state. In contrast to the case with- 
out driving, the manifold has bounded effective barriers 
or a finite typical barrier at finite velocity (maximum ef- 
fective barrier). In a stationary state, where the manifold 
is driven by a force F, the typical barrier height U deter- 
mines the rate of jumps in direction parallel or opposite 
to the driving force. The resulting velocity v is estimated 
by 
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-U/T 



liF (7) 



for F£_ < T with an effective mobility /x = £, 2 /t T. The 
scale £ is related to the size of the barrier and therefore 
depends implicitly on U. However, this dependence is 
much weaker than the exponential dependence displayed 
explicitly in Eq. (7) and can be ignored. 

Inverting Eq. (7) one can define from the stationary 
velocity-force characteristics (VFC) v s t(F) a typical bar- 
rier for stationary creep as a function of the driving force 
by 



U Bt (F)=T]n\ f iF/v et (F)]. 



(8) 



The scaling theory predicts a dependence of the bar- 
rier height on the driving force according to U s t(F) ~ 
F -e/(2-0 g. 

The structure of the manifold in the driven state can 
be determined as follows: the typical barriers are the 
maximum effective barriers, which are able to pin local 
regions of the manifold temporarily. The manifold will 
be rough according to Eq. (g) only up to the size of these 
regions. This size is related to the maximum effective 
barrier U s t via Eq. (||). 

If the structure of the manifold initially does not cor- 
respond to the stationary state, it will relax with a re- 
laxation time r(U si ) given by Eq. (||), since U st is the 
maximum effective barrier. However, this time scale is 
relevant only very close to the stationary state. 

Heuristic model. We now develop a heuristic model 
describing the manifold dynamics far away from the sta- 
tionary state and the response of the manifold to a time- 
dependent driving force. In this case the manifold needs 
a finite time to adapt its structure to the varying force. 
For simplicity we describe the structure of the manifold 
by a single parameter, the length scale L c (t) up to which 
barriers are effective, or equivalently the corresponding 
energy scale U(t). We assume that the velocity responds 
instantaneously to the force according to 



v{t)= i iF{t)e- u ^' T . 



(9) 



This means that the structure at the time t determines 
only the effective mobility at that moment. 

Due to the finite relaxation time, which grows infinitely 
for the decreasing velocity, the structure itself will not 
adapt itself instantaneously to the driving force, but it 
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will rather depend on the history of the system. We de- 
scribe this history dependence effect by a simple equation 
for the evolution of the typical activation barrier U(t): 



dt w t 



-U. t (F(t))/T 



(10) 



This is the central equation for our study of hysteretic 
creep. The first term on the right hand side represents 
the tendency of the manifold to get pinned by the dis- 
order potential. Thereby it gains potential energy and 
and the barriers increase. This tendency is balanced by 
the second term, which incorporates a saturation of the 
typical barrier due to thermally activated depinning of 
local regions for finite force. 

Although a microscopic derivation of Eq. ([n]) is be- 
yond our capabilities, we would like to argue that it is 
the simplest conceivable equation of motion of U which 
satisfies the following conditions: (i) in the absence of dis- 
order one has an ohmic stationary VFC with U s t (F) . The 
equation of motion is then solved by U(t) = indepen- 
dent of the driving force, (ii) In the absence of a driving 
force, U(t) must grow logarithmically as in Eq. (||). This 
is indeed the case, since e _c/st ' i?=0 ^ T vanishes in a glassy 
state, (iii) When U(t) — U a t(F(t)), the manifold is in in- 
stantaneous equilibrium, where U (t) must not change, 
(iv) If the driving force F > is constant in time, the 
structure has to approach its stationary state with a fi- 
nite relaxation ti me (|3] ). This behavior is readily verified 
by linearizing Eq. (|1C|) with respect to U(t) — U st (F) . The 
Eq. (10) may be also considered as a simplified version of 
the evolution equation for a barrier distribution function 
in glassy systems used in Q. 

Phenomenology of the model. In general, the hysteretic 
behavior emerges as the system parameters (like the driv- 
ing force and/or temperature) change on a time scale of 
the order of or smaller than the characteristic relaxation 
time. Hence we expect these hysteretic phenomena to 
be particularly pronounced in glassy systems, where the 
equilibrium relaxation time diverges and aging effects oc- 
cur. However, since in the driven case the external force 
F sets the upper limit U s t (F) for the barriers the system 
can be trapped into, the relaxation times are large and 
we expect noticeable effects only at sufficiently small ve- 
locities. As we shortly show the equation of motion jl0| ) 
is not only able to capture hysteretic effects but even 
more it exhibits a retarded and almost threshold-like on- 
set of motion, when a driving force is switched on after 
the system has been aged at zero driving force. 

We performed a numerical integration of this differen- 
tial equation for some time-cycles of the driving force. 
These calculations have been performed for a station- 
ary characteristics with U s t(f) = T/~ - 5 . Here / is a 
dimcnsionless force given by the ratio of F and the de- 
pinning force. The scaling exponent is that for the three- 
dimensional vortex lattice with logarithmic roughness, 
C = |1C[]. Time is measured in units of to- 
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FIG. 1. Depinning of the manifold, when a constant driv- 
ing force is switched on after aging at zero drive. Inset: Time 
evolution of the typical barrier U(t). Main plot: Time evolu- 
tion of velocity v(t). The curves (from the left to the right) 
correspond to initial barriers U(0) =125, 250, 500, and 1000. 

First we examine the time dependence of velocity in a 
depinning scenario. It is supposed that the manifold has 
initially a large typical barrier U(t = 0), which can be ob- 
tained by aging at zero drive. Then at t = a finite drive 
/ > is switched on such that U st (f) < U(0). This leads 
to a very abrupt onset of motion illustrated in Fig. [l]. 
This behavior can be understood as follows: initially the 
barrier decreases linearly in time. This implies that the 
velocity increases exponentially in time. The time scale 
is set by the relaxation time of barriers of height U s t (F) . 
After a time t dcpin rj T(U st (F))U(t = 0)/T the barrier 
height and drift velocity reach their stationary values. 
On a linear velocity scale the exponential increase and 
the saturation appear as an almost discontinuous onset 
of motion. 

The sharpness of the onset of motion increases for 
smaller forces, since there the stationary VFC has a 
stronger nonlinearity. In the example of Fig. |l| the re- 
laxation time r(C/ s t(/)) ~ 2 • 10 5 is smaller than tdcpin by 
a factor of order 10 2 . The depinning time tdcpin, which 
is essentially proportional to the initial barrier, increases 
only logarithmically as a function of the aging time, dur- 
ing which the manifold has to be kept at zero drive in 
order to increase the initial barrier height. 

Next we perform sweeps of the driving force in time in 
order to demonstrate the possibility of an extreme hys- 
teresis in v(f). As a generic example we consider the 
case, where an initial force (/ max = 6 • 10~ 3 ) well below 
the depinning force (/ = 1) is switched off with a con- 
stant rate (during a time t SW ccp = 10 10 ). The sweep time 
is chosen to satisfy r({7 st (/ max )) < t 
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The inset of 

Fig. H illustrates the time development of U (t) during 
this sweep. It increases first linearly and then logarith- 
mically, since finally the system ages as in the absence 
of a driving force. Then the force is ramped up to / max 
again with the same rate, immediately after / = has 
been reached. The barrier continues to grow until the 
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driving force becomes so big that U st (f(t)) > U(t). In a 
linear plot as in the main part of Fig. 0, the characteris- 
tics exhibits practically no hysteresis, i.e. v(t) coincides 
with u 8 t (/(*))• 

Hysteresis appears only after the system has been aged 
at zero (or very small) driving force. In the numerical ex- 
ample depicted in Fig. || we suppose that the system has 
been aged a certain time until the barrier reaches values 
U = 125, 250, 500, or 1000 (full lines from the left to the 
right). These values correspond to stationary barriers at 
forces of order / ~ 10~ 6 <C /max- Starting from these 
barrier heights at t = 0, the force is switched on up to 
/max linearly during i SWCC p = 10 10 . Initially the barri- 
ers decrease extremely slowly and then faster and faster, 
until U(t) w C/ S t (/(*))■ This again leads to a threshold- 
like onset of motion from an effectively pinned state to 
a state, where the manifold moves with the velocity v(t) 
located on the stationary creep VFC. As we increase bar- 
riers U(t) reflecting the growing waiting time at / = 0, 
the threshold force /th(i) shifts upwards to larger values. 
Note the striking similarity of the obtained hysteretic 
threshold-like onset of motion to the hysteretic switching 
that has been observed in recent experiments on high- 
temperature superconductors M. 

In conclusion, we have presented a simple phenomeno- 
logical model of elastic manifold driven through quenched 
disorder which displays hysteretic dynamic behavior. Dy- 
namic freezing of the system into deep metastable states 
at zero driving force leads to a threshold- like onset of mo- 
tion at the threshold force f(t w ) growing with the wait- 
ing time t w at / = 0. It is important to note that the 
presented threshold dynamics is different from the usual 
depinning transition at zero temperature. The threshold- 
like onset of motion demonstrated here is specific to finite 
temperatures and occurs well below the depinning tran- 
sition. It would be interesting to perform a systematic 
study of the threshold dynamics in the vicinity of the 
melting transition, where this dynamic hysteresis inter- 
feres with the static hysteresis of melting proposed in || . 
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FIG. 2. The driving force is swept down from / = 6 • 10 -3 
to / = with constant rate in time t = 10 10 . Then it 
is immediately swept up again with the same rate. Ini- 
tially U = U s t{f) is supposed. Inset: U(t) during sweep 
down for < t < 10 10 and during subsequent sweep up for 
10 10 < t < 2 ■ 10 10 . Main plot: v(f) shows no noticeable 
hysteresis. 
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FIG. 3. The driving force is swept up from / = to 
/ = 6 ■ 10~ 3 with constant rate in a time t = 10 10 . Different 
initial conditions are used, see text. Inset: Time evolution of 
U(t). The larger U initially, the later it comes down to the 
stationary value. Main plot: v(t) versus f(t). The larger U 
initially, the later v jumps up to the stationary VFC (dashed 
line) . 
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